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Abstract 

A theory of integer quantum Hall effect (QHE) in realistic systems based on von Neu- 
mann lattice is presented. We show that the momentum representation is quite useful and 
that the quantum Hall regime(QHR), which is defined by the propagator in the momentum 
representation, is realized. In QHR, the Hall conductance is given by a topological invariant 
of the momentum space and is quantized exactly. The edge states do not modify the value 
and topological property of a xy in QHR. We next compute distribution of current based on 
effective action and find a finite amount of current in the bulk and the edge, generally. Due 
to the Hall electric field in the bulk, breakdown of the QHE occurs. The critical electric field 
of the breakdown is proportional to I? 3 / 2 and the proportional constant has no dependence 
on Landau levels in our theory, in agreement with the recent experiments. 



1. Introduction 



We study the quantum Hall effect in realistic systems in the present paper. Especially 
finite size effects and finite current effects are investigated in details. 

Two-dimensional electrons in a strong perpendicular magnetic field have discrete en- 
ergies with a magnetic field dependent finite degeneracy per area at each energy. When 
degenerate electrons are distributed uniformly distance between two electrons has a mini- 
mum value. Conjugate momentum, hence, is defined on a torus. Thus the momentum space 
is compact. The energy space seems not to be compact but physical requirement to the 
propagator at po = ±00 or ±200 leads the energy space defined by the propagator compact. 
A representation of having these properties in a manifest manner was constructed based 
on von Neumann lattice of guiding center coordinates, and exact low energy theorems were 
given based on this representation 1,2 . Namely, it was shown that Hall conductance a xy 
is a topological invariant of a mapping from the momentum space to a space defined by 
the propagator, and the quantization of the a xy as (e 2 /h)N at the plateau, thus, is exact 
in systems of interactions and disorders. One-particle properties in systems of short range 
impurities, boundary potential, and periodic array of potentials have been obtained 3 also 
based on the same representation. Localization is shown easily in our method. 

From our previous investigations, the a xy becomes a topological invariant of the mo- 
mentum space and is quantized as (e 2 /h)N exactly at the plateau if all the one-particle 
states around Fermi energy are either localized or have finite energy gap. Hereafter we call 
this energy region as quantum Hall regime(QHR). In finite system with boundary, there are 
edge states 4 which are extended along the edge and have continuous energy across the Fermi 
energy, generally. They may give a finite correction to the quantized value. We will show 
in the present work that a resummation formula of the a xy with the momentum represen- 
tation is valid in finite systems as well and QHR is defined based on the propagator in the 
momentum representation. The outside regions of the continuous band of propagator in the 
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resummation formula agrees with QHR. Perturbative series converges well in QHR. Conse- 
quently, a finite size correctoin disappears in QHR. Niu and Thouless 5 claimed before that 
the finite size correction is of exponential type. Our results show that even an exponential 
correction does not exist in QHR. 

The quantized Hall conductance 6 is used as a standard of resistance and to determine 
fine structure constant, a. Its precise value is needed in testing quantum electorodynamics. 
The theoretical foundation of the quantum Hall effect (QHE) is, however, still insufficient. 
Especially there is a controversy in current distribution in a finite system and in finite size 
corrections of the QHE. Breakdown of the QHE is also observed recently if the current 
becomes large. 

Experiments are done with semi-conductors of Hall bar geometry. There are edges in 
such systems. Concerning edge states, a controversy is the following. In one method for a 
proof of QHE, Biittiker-Landauer 7 formula is used. It is assumed that the one-dimensional 
edge states are the only current carrying states around the Fermi energy and are connected 
with the leads. The quantization of the Hall conductance due to the one-dimensional edge 
states was derived. A correction may be caused by backscattering of edge states. In other 
approaches, the bulk states carry the currents and their roles are impotant. The edge states 
might give a correction to the quantized value, because they are extended in one direction 
and have continuous energies. The role of edge states and bulk states are thus reversed. 
Two approaches are thus very different. It is a purpose of the present work to resolve these 
controversies and to give a foundation of the quantum Hall effect in finite systems. 

We study a finite size effect and a finite current effect based on von Neumann lattice 
representation in the present paper. We show that the momentum representation and re- 
summations of the diagrams based on the momentum representations are valid in finite 
systems and the QHR is defined as the outside region of the continuous energy band using 
the resummation formula of the a xy . It is shown that the quantization is exact in QHR of 
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finite systems despite the fact that the edge states have continuous energies and cross Fermi 
energy. There are current flow both in the bulk and near the edges generally and due to the 
Hall electric field at the bulk, the QHR becomes narrower as the current becomes larger and 
eventually disappears at a critical value. QHE is broken, then. We estimate a critical Hall 
electric field and find an agreement with the experiments. 

The paper is organized in the following manner. We review our representation of two- 
dimensional electrons in the magnetic field based on magnetic von Neumann lattice in the 
rest of Section 1. In Section 2, one-particle properties are studied. Exact low energy theorem 
concerning the slope of current correlation function in infinite system are given in Section 
3. Topological invariant expression of Hall conductance is given and the corrections due 
to interactions and disorders are shown to disappear in QHR. Finite system with cylinder 
geometry is discussed in Section 4 and finite systems with Hall bar geometry, which is the 
geometry of realistic experiments, are discussed in Section 5. In Section 6, the current dis- 
tribution is computed with a use of effective potential. Breakdown of QHE is also discussed. 
A critical value of Hall electric field from our theory is proportional to _B 3 / 2 and the propor- 
tional constant is independent of Landau levels. These properties are in agreement with the 
recent experiments by Kawaji et al 8 . Summary is given in Section 7. 

We review our representation 1 ' 2 ' 3 of two-dimensional electrons in a strong perpendicu- 
lar magnetic field here. It is based on von Neumann magnetic lattice 9 and has excellent 
properties for a purpose of studying one-particle properties and of giving the exact low en- 
ergy theorem of quantum field theory in quantum Hall system. Localized states are studied 
with the coordinate representation and extended states are studied with the momentum 
representation. 

A one-body Hamiltonian of a planar charged particle with a strong perpendicular mag- 
netic field, 
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H = M±^!, dlA2 - d 2 A 1 = B, 
2m 



(1.1) 



is expressed as 



* = ^«W), (1-2) 



with two-dimensional relative coordinates in the magnetic field 10 . They are defined by 



f = ^(Py + eA y), 

V = —^{Px + eA x ), (1.3) 



and commute with another set of coordinates, center coordinates, (X, Y) defined by 

x = x -i, 

Y = y-r), 

ih ( L4 ) 

K,JC] = K,y] = [i7,x] = [i7,y] = o. 

It is convenient to use variables (X, Y) and (£, 77), instead of (x, y) for studying the electrons 
in the magnetic field. A set of well localized functions, 



fl(C,v) ® |-Rm,n), 

(1.5) 

/ifo^) ® (-Rm.nl, 



defined by 



m 1 



|jj> ^ _ ^_^mn+m+n e A 1f y/ir(m+in)— A^/iv(m—in) |q^ 

I jd n \ _ „i7r[(mi-m2 + l)(ni-n 2 + l)-l] „-(7r/2)[(mi-m 2 ) 2 +(ni-n 2 ) 2 ] 

\ Jl mi,ni |- Il m2,n2/ — c c 5 

(Rm 1 ,n 1 \Rm2,n 2 ) = (^mi,m |-^-m2,n 2 ) = 0; ^ 



mi, ni m 2 ,n 2 



«m,n = a(m,n), m,n : integer, a = W — — , 

A|0) = 0, A = ^(X + tY), 
A\ R m,n) = ^(m + in)\R m , n ), 



is a complete set and is used as base functions. Dual basis (i2 m , n | is defined, with a Green's 
function G(mi, n\\ m,2, 77,2), by 

|= G(mi,ni;m2,n 2 ){R m 2 ,n 2 1 5 

m 2 ,n 2 



m,n 

G(mi, ni; m2, 712) is well localized in small |mi — 7712 1 and |ni — 722] region as is shown in 
Ref.(3). We expand the electron field as, 



* = Yl a l( m > n )M^V) ® |-Rm,n), 

'' m ' n (1.8) 

m,n\i 

l,m,n 



and regard coefficients ai(m,n) and bi(m,n) as quantized operators. The free kinetic term, 
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potential term, and the electromagnetic current are expressed with these new variables as 

ix^\x) ty(x) = Eib(m,n)a(m,n), 

2m 

l,m,n 

J dx^(x)^(x)V(x) = ^ b h {mi,ni)ai 2 (m 2 , n 2 )Vi 1 ,i 2 (mi, n\\ m 2 , n 2 ), 
J n( x ) = ^ bujmi, ni)a; 2 (m 2 , n 2 )r^ 2 (mi, m; m 2 , n% x), 

Vi u i 2 ( mi ,ni;m2,n 2 ) = j ^^{R mi ,n,\e ikX \R m2 ,n 2 )Ui^ iki %)V(k), ( L9 ) 

T l ^(m 1 ,n 1 ;m 2 ,n 2] x) = J ^(^.nJc^li^^J^^e^+^/iJe-^, 

^ = (1, -{eB/m)ri, (eB/m)S), 
Advantages of the present representation are summarized in the following: 

(1) The kinetic term is diagonal in the Landau level index and two-dimensional lattice 
coordinates and only the potential term has off diagonal terms. Hence the localization 
by short range random impurities in energy regions \E — E\\ > 5 with a small impurity 
dependent parameter 5 are shown easily in the present method. 

(2) Functions used as basis and dual basis are well localized around center coordinate 
Rm t n. Hence the multi-pole expansion of current operator itself and of commutation relations 
between the current and field operators are applicable. Due to these properties it is a 
straightforward matter to derive Ward-Takahashi identity and related identities as well as 
the exact low energy theorem. Coordinate dependence of a potential V(x) is preserved in 
Vj^^mi, m 2 ; n\, n 2 ) 9 with a smeared form. For instance, a short range potential, 

V(x)=g5(x), (1.10) 

is transformed into 

Vb, (m 1 ,m 2 ;n 1 ,n 2 ) = 4e" f(m?+ri?+Jni+ni) [l + O(e-i)], (1.11) 

CL 

which is a short range potential with an extention of few magnetic distance. Other compo- 
nents Vi 1 j 2 (mi,m 2 ;ni,n 2 ) have the same properties. 



(3) Apart from disorder potentials, translational invariance is manifest and Fourier trans- 
formation of the operators can be denned, because magnetic translations commute each other 
and constitute abelian group in our lattice system. Extended states can be studied with the 
momentum representation. 

(4) The coordinates are defined on lattice sites. Since the momentum defined by Fourier 
transformation from lattice sites is defined on a torus, which is compact, the a xy becomes a 
topological invariant in momentum space manifestly. Neither artificial boundary condition 
in configuration space nor periodic potential are needed. 

§2. One-particle states in finite systems with disorders and interactions 

We study one-particle properties of systems with one short range impurity, dilute short 
range impurities, boundary potential, and interactions between electrons in this section. 
Localized states are studied with the coordinate representation and the extended states are 
studied with the momentum representation. It is shown that the QHR which is defined by 
the propagator in the momentum representation exists in finite system if the magnetic field 
is of suitable strength. 

(2- a) Disorder potentials 

A Hamiltonian, 

J rf^t( g )[ (P+^) + V (x)Mx) = b h(Ri)[E h 5 h ,lJm,R2 + V ll!h (R l ,R 2 )}a l2 (R 2 ) 

h,h,Ri,R2 

(2.1) 

discribes the systems with a disorder potential V(x). The kinetic term becomes diagonal 
and the potential term becomes non-diagonal. In fact this term becomes a lattice kinetic 
term if a suitable periodic potential is used for V(x), as was shown in Ref(3). Eigenstates 
are extended, then. For a short range potential, on the other hand, eigenstates of energy E 
in region, \E — Ei\ > 5 with a small parameter 5, becomes localized. This is easily seen from 
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the eigenvalue equation, 



£ V lul2 (Ri,R2K\R 2 ) = - ^>n£i). (2.2) 

Eigenvector u^\mi,n\) has the same coordinates dependence as the transformed potential 
term V^j^Ri, R2) provided E — E\ 7^ 0. Hence an eigenvector is localized around a short 
range impurity if its energy E is away from Landau level energy E\. We confirmed this 
property of wave functions by solving the equation numerically 3 . Wave functions of the 
energy E — E\ 7^ have spatial extensions of few magnetic lengths. 

Many short range impurity problem is studied as easily as one impurity problem if 
impurities are dilute and random. Localized wave functions around one short range impurity 
decrease very fast and almost vanishes at nearby impurities if distances between impurities 
are much larger than the magnetic distance. Corrections of localized wave functions due 
to nearby impurities are thus very small and is treated perturbatively. If E — E\ is also 
very small, the situation is different and the energy denominator becomes very small and 
compatible to the small numerator in perturbative series. The correction can become order 
one, and perturbative series may be divergent, then, and wave function becomes extended. 
Consequently, one-particle property depends on its energy in the system with dilute random 
impurities. Wave function is localized if its energy E is in \E — E\\ > 5 with impurity 
dependent small parameter 5. The propagator in the momentum representation has no 
singularity in this energy region. 

(2-b) Boundary potentials 

A potential Vq6{x) shows a potential barrier in the positive x region with a boundary 
at x = 0. Due to a translational invariance in ^/-direction, it is convenient to project the 
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potential and wave functions to those of definite momenta, as 



a l" K ' a 

V hM {m 1 ,n 1 -m 2 ,n 2 ) = — / dpye ip y {ni - n2)a V llh (m 1 ;m 2 ;Py), 

^ J-n/a 

ui (m, n) = ^ — J dp y e tp y na ui (m, p y ). 



(2.3) 



The eigenvalue equation becomes, 

2^mi,m2 + Vi lh (m 1 ;m 2 ;py)}u < lf(m 2 ,py) = E^u\f(mi,p y ). (2.4) 

h,m 2 

We solved the equation numerically. Edge states have continuous energy and wave functions 
are extended in y-direction and are confined in x-direction. As the boundary potential gives 
steep electric field in x-direction, edge states have current with definite direction and are 
regarded as chiral modes. They have been computed in the present representation and 
were given in Figs. 7 and 8 of Ref.(3). They have continuous energies in the wide range and 
bridge from one Landau level to next Landau level. Hence there are always one-dimensionaly 
extended states around an arbitrary value of Fermi energy. 

In systems with finite width in the x-direction, there are two edges. The wave functions 
change their properties depending upon their widths. To see them, we solve systems of 
potential wells in the x-direction numerically. As are shown in Fig.l and Fig.2, if the width 
is large, two edge states are separated, but if the width is small, they are combined. In 
wide systems, one-dimensionaly extended states have small overlapp with the momentum 
eigenstates that are used in our method of computing the Hall conductance. We will see 
that they give no correction to the quantized Hall conductance in certain situation. 

Now we find out when the edge terms and impurities give no correction. In the section 
5 of the present paper, we apply the momentum representation of the a xy as in Ref.(l). The 
edge potential term is treated as perturbative term in the Hamiltonian and in the propagator. 
In this expansion, infra-red divergence is not involved unless the Fermi energy agrees to the 
Landau level energy. Moreover, the series converges well if the energy difference, E — Ei, 
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is much larger than the perturbative energy of the states projected into the momentum 
eigenstates. These momentum eigenstates, although they are fictitious, play the important 
roles in our calculations and a band width of these fictitious extended states is important and 
is estimated in the following. If the Fermi energy is in the outside of this fictitious band, the 
convergence of the perturbative expansion is very good, and the low-energy theorem for the 
<j xy is derived easily. It should be noted that this can occur even if the real one-dimensional 
edge states bridge from one Landau level to next Landau level and have zero energy around 
the Fermi energy. 

For estimation of the band width in the momentum representation we study eigenvalue 
equation of the Hamiltonian in one method, and we study the propagator in another method. 

The band width in the first method is estimated directly from Eq.(2.2). The boundary 
potential in magnetic lattice representation V^ ^Ri, R2) decreases rapidly with the distance 
between the coordinates R\ or R\ and the boundary and behave as 

V(Ri, R 2 ) ~ VoMax(e-^) 2 , e~i ( ^ )2 ), (2.5) 

where Vq is the potential height and L\(L2) is the distance. Hence the equation(2.2) becomes 

— * 

at a point R in the middle of the system, 

.(«)/ 



^1 I ( . — *. . ( ru\ . — * 



\u£\R)\ \ul a \R)\ 



(2.6) 



An extended eigenstate which does not vanishes in whole space satisfies, 

\u^\R)\ > -fL (2.7) 
By combining Eq.(2.6) and Eq.(2.7), we have 

-Ei\< cV VLe-^^ )2 . (2.8) 
The band width, A e d ge , thus satisfies, 

Acdgc < cVoVZe-^^ 2 . (2.9) 
The extended states in the momentum representation are generated from the extended states 
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which have energies within Eq.(2.8). 

In the second method, we study the propagator. As is expressed in Ref.(l), the energy 
of the momentum eigenstates in the system of disorders can be defined and is expressed 
from diagrams of Fig.3. We calculate the lowest non-trivial order contribution here. The 
characteristic properties of the full order are known from them. The self-energy correction 
of the momentum p up to 0(V 2 ) is given by 



where the state has definite momentum and the boundary potential of Eq.(5.5) and impurity 
potential are used. Since the boundary potential is not vanishing only in the finite region 
around the boundary of few magnetic length, the above quantity has the magnitude, 

^ + ^M^) 2 , (2.11) 
L E — Ei L ' V ; 

with magnetic field independent constant c and a width L. If the magnetic field is strong 
enough, or the width L is large enough, then the self-energy correction due to the boundary 
potential or due to the boundary potential and the impurity potential becomes arbitrary 
small. Hence the width of the fictitious band which appears in our computation can be 
(much) smaller than the Landau level spacing if the L is large enough or the magnetic field 
is strong enough. If the Fermi energy is in the outside of this fictitious band, the system is 
in the QHR. The QHR can be realized in the finite systems. 

(2-c) Electron interactions 

Interactions also remove the degeneracy of Landau levels. We estimate band width due 
to Coulomb interactions based on perturbative calculation for momentum eigenstates of 
Fig. 4. Within the Z-th Landau level space, the propagator is written as, 



where v is the filling factor in the Z-th Landau level. The eigenvalue of the energy is obtained 
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(2.10) 



1 



po- Ei- E(p,po)' 



(2.12) 



m="J J^ 2 V(k)\(fl\e ik W, 



from the pole of the above propagator. The energy shift from the degenerate value E\ is 
given by £(p), which vanishes at v — 0, and agrees with 

AE = I J^ V ^\ eik W = z 2 ^B- (2-13) 
at v — 1 in the lowest Landau level space. The difference of the energy correction between 
v — 1 and v — gives the band width due to interactions. We have thus 

= eV27re£, (2.14) 

which has a weaker magnetic field dependence than Landau levels spacing. 

The width in our momentum representation due to edge states, Eq.(2.9) or Eq.(2.11), 
becomes (much) smaller than the Landau level spacing with a strong magnetic field or with 
a large L x . The width due to interactions, Eq.(2.14), also becomes much smaller than the 
level's spacing. Hence with a strong magnetic field, there are wide energy region which have 
no singularity due to two-dimensional extended states. The localized states due to short 
range impurities can have energies in these regions. We call these regions as quantum Hall 
regime(QHR). QHR exists in systems of boundary, disorder and interactions if the magnetic 
field is strong enough. Whole spectrum is shown in Fig. 5. 



]3. Low energy theorem of a 



xy 



Low energy theorem of a xy is satisfied in QHR with impurities and interactions. We 
study many-body problems with (quantized) field operators {a/(m, n; t), bi(m, n; t)} which 
are defined on lattice sites, 



2nh 

R mn = a{m, n), a = d — (3.1) 
and satisfy an equal time commutation relation, 

{ai 1 ( y m 1 ,n 1 ;t 1 ),bi 2 (m2,n 2 ;t 2 )}5(t 1 -t 2 ) = Si u i 2 5m 1;rn2 5 ni! n 2 5(ti -t 2 ). (3.2) 
Operator of zero momentum was not included originally in operator set, but is added for 
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computational convenience 1 . This operator decouples from physical space due to the con- 
straints of von Neumann lattice coherent states. 

The propagator, vertex function, and current correlation functions in translationaly in- 
variant systems are defined by, 

/ dhdt 2 e i(pi " i ^ 2) (n«/ 1 (^i^i)^(^2,t 2 )}) = 4( 27r ) 35 (w-P2)^ 1 ,/ 2 (w), 

/ CL 

Xi ,X2 

I dt x dt 2 dx e l ^ lXl ^ 2X2 )(T{j M (a:)a; 1 (X 1 ,t 1 )^(X 2 ,t 2 )}) = 

J x u x 2 (3.3) 



a 



-(2n) 3 5( Pl -p2 + q)Si lA (p 1 )T l ^(p h p 2 )S l y 2 (p2), 
J dx 1 dx 2 e i ^-^(T{j fI1 (x 1 ) 3fl2 (x 2 )}} = (2n) 3 5(q 1 ~q 2 )ir lll ^( qi ). 

T(-) denotes time-ordered product and momenta conjugate to lattice coordinates are defined 
on a torus. Hall conductance in QHR is computed from the time ordered current correlation 
functions and agree to that of retarded product used in Kubo formula 11 . Owing to the 
current conservation, 

d„f = 0, (3.4) 

and the commutation relations Eq.(3.2) and their representation Eq.(1.9), the above Green's 
functions satisfy 1 ' 2 , 

Q^u(q) = ^iiAq)q v = 0, 
<f r>i,p 2 ) = S-\ P1 )R( P2 ) - L( Pl )S~\ P2 ), 

(3.5) 

Rh,iM = 5 hh + nx[d x {p)6 hh + d xhh ] + iq y [d y (p)5 hM + dy M2 ], 
L h ,h(p) = 5 hh + iqx[d' x (p) S h,h + d' xhh ] + iq y [d'y(p)6i lM + d' yMM ], 

di {p) = dj(p) = e ip(Xl ^ 2) (^i|(^-^i)^2), 

Xi-x 2 (3.6) 

dihhip) = d' ihh {p) = (fh\(i)i\fi 2 )- 
The explicit forms of di(p) and di(p) are given in Appendix A. The Hall conductance is the 
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slope of n flu (q) at the origin and is written as 



>xy 



d 



3! dq p 



q=0 



2ir 



N w , 



W 247T 2 

with the transformed propagator 



1 f, 3 Ul/prjl l dS' 1 (p) d , ,dS- 1 (p) d , ,dS- 1 (p) d , 



(3.7) 



dp p 



dPv 



dp v 



S(p) = V(f)S(p)U(p), 

U(p)^U- 1 (p) + l i (p) = 0, 
opi 



dV~ l {p) 
dpi 



(3.8) 



V{p)+r i { P ) = 0, 



kip) = ri{p) = di(p)8 hh + di{p)8 hh . 
We have used a fact that the vertex function which is transformed by the above matrices as 



f>i,p 2 ) = u- 1 (p 1 )r fl (p h p 2 )v~\p 2 ) 



satisfies the standard Ward-Takahashi identity, 



f M (p,p) = 



dp P 



(3.9) 



(3.10) 



It should be noted that the relations Eqs.(3.4)~(3.10) are satisfied in interacting systems as 
well. 



The meaning of N w : 

The N w of Eq.(3.7) is a winding number of mapping from the momentum space to a 
matrix space defined by the propagator Si 1; i 2 (p). The momentum in spatial direction is 
defined on a torus. The integration region of the momentum in temporal direction is also 
regarded as a compact space, when the propagator is non-singular at po = ±oo. Physically, 
there is no singularity at po = ±oo in theories we study here and magnetic field effect should 
disappear at po = ±oo. Furthermore, the energy integration region becomes a closed path in 
complex plane, which is compact, if the difference N w (Ep 1 ) — N w (Ep 2 ) is computed. Hence 
the momentum in temporal direction is also regarded as compact. 
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Consequently, the momentum space is regarded as a compact space. The matrix space 
is of infinite dimensions and seems to have a difficulty in defining a topological invariant. 
However, since the energy of Landau level E\ diverges if I — > oo, matrix space which corre- 
sponds to finite energy is of finite dimensional. It includes SU(2) space as a subspace. Hence 
N w generally agrees to an integer as an elements of 113(G), G D SU(2). 

Computation of N w '■ 
(i) Free system 

We compute N w in a free system without disorders and interactions first. Let substitute 
Eq.(3.6), (3.8), and free propagator Sq(p) expressed by (po — ±ie)~ 1 5i 1 j 2 into Eq.(3.7), 
then we find that the integrand of Eq.(3.7) is independent of p and is given by 



N w depends on the Fermi energy and agrees to an integer at E 7^ E\ as is shown in Fig.6. The 
value is ambiguous at E = E\ due to degeneracy of Landau levels. By adding momentum 
dependent small energy, the degeneracy is removed and we have a unique value of N w . The 
value thus obtained is proportional to electron filling factor u, as 




(3.11) 



We have, then, 



N w = l, Ei<E F <E l+1 . 



(3.12) 



(3.13) 



which leads the following classical value of the Hall conductance: 



a 



xy 




(3.14) 



(ii) Systems with impurities 
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Degeneracy of Landau levels are partially removed by impurities and there appear one- 
particle energy levels of having energy E at E ^ E\. It was shown that the energy eigenstates 
of E at \E — E\\ > 5, are localized if the impurities are dilute, random and short range based 
on perturbative expansion method in von Neumann magnetic lattice representation, where 
5 is a system dependent small constant. Numerical calculations in random system also show 
that the wave functions are localized at E ^ E\ with the localization length being inversely 
proportional to some power of E — E[, (E — E[)~ x , A > 0. Since the wave functions 
of localized states have finite extensions and their energies are discrete, it is possible to 
compute the contributions of the localized states to a xy perturbatively, if the Fermi energy 
is in the range \Ep — E\\ > 5. 

We assume, hereafter, that the system is in quantum Hall regime where all the one- 
particle states around Fermi energy are localized and the momentum conserving propagator 
in the momentum representation has no singularity. We show that only the momentum 
conserving amplitude contributes to the linear coefficient of the current correlation function 
first. We decompose the current correlation function into the momentum conserving part 
it${q) and the momentum non-conserving part {qi, Q2) ■ 

o (2) 

n>Wi>92) = (2?r) 5(qi -92)^(91) + ^U(qi,Q2)- (3.15) 

In the latter, q\ and q\ are momenta carried by currents and are independent each other. 
This term is not zero generally due to a lack of translational invariance. Current is conserved, 
hence the identities, 

(3.16) 

9if $(91,92) = ^$(91,92)92 = 0, 

are satisfied. In quantum Hall regime, n^u(o) and tt^}(qi,Q2} have no singularity at the 
origin = or q\^ = q^ = 0, because one particle states have discrete energies and their 
wave functions have finite spatial extensions. Singularity is not generated. Thus they are 
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expanded as, 

*$(q) = 7r#(0) + Trgp(O)^ + 0(g 2 ), 

(3.17) 

g 2 ) = 45(0, 0) + 7rg p (0, 0)gf + tt^Uo, 0)g£ + 0(g 2 ). 
The coefficients are finite. We plug these forms to Eq.(3.16), and we have, 

45(0) = 

vtSp(0) + 4pV(0) = 0, (3-18) 
4i;p(0) + 4^(0)=0. 

The slope of vr^(g) at the origin is, hence, totally anti-symmetric, and is written with one 
constant c, as 

We have, also, 

45(0,0) = 0, 

4?;p(0,0) = 0, (3.20) 
4^(0, 0) = 0, 

by substituting the second form, Eq.(3.17), into the second equation of Eq.(3.16), because 
qi and q\ are independent. Thus the linear term vanishes if q\ and q\ are different and 
independent. In fact, even in systems with translational invariance in one direction where 
one component of qi are the same, the linear term vanishes. This is understandable easily 
from the fact that there is only one free parameter left in Eq.(3.16), and one additional 
condition makes non-trivial solution disappear. The slope of the current-current correlation 
function at the origin is due to the momentum conserving part, Tr$(q), and the momentum 

(2) 

non-conserving part, 7v^J(qi, ^2), does not contribute. This will be used also in the next part 
when interaction effects are studied. 

(iii) Resummation formula of the a xy in the momentum representation. 
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We use a representation of the momentum conserving part, nj^j(q), in terms of a mo- 
mentum conserving part of propagator and a momentum conserving part of vertex part. 
They are defined by summing the momentum conserving parts in perturbative expansions. 
The details of the resummation formula have been given in Ref.(l). Here, we give only the 
momentum conserving part of propagator in a system of impurity potentials or boundary 
potentials as an example. QHR is defined as the outside region of the energy band of the 
momentum conserving part of the propagator. In order to write the propagator with a 
momentum representation in the system of disorder potential V, we start from, 

, | 1 | 1 rr (Pl|V>2) (gljjV) V\V\P2) , 

E — (H + V) m) E - £b(pi) [ PUP2 E- £b(pi) E- E ( Pl ) E - E (p>) + " ' J ' 

(3.21) 

H \p) = E (p)\p), 

where Ho is assumed to be invariant under translations. Eigenvalue Eo(p) could have a 
dependence on the momentum. In quantum Hall system, Hq is given by 

H = Y t Eibl(R)a l (R), (3.22) 

and the energy eigenvalue, Eq(p), thus is constant, E\. Due to the non-invariant term, V, the 
momentum is not conserved and p2 is different from p\ generally. The momentum conserving 
propagator within the Z-th Landau level space is defined by the perturbative series as, 

= (PI^^IP) = ^[1 + + Te (3^3) 

Hi — ^X2 Q + V ) Hi — hil *-> — — till tli — Hil 

This can be written in the following form: 

s{p) _ 1 {1 , (P\V\P) , y WW) (P'\V\P) , 

p>=p 



| y {p\V\ P ') {p'\V\ P ) | ] 

p'^p 

TP Z?L _,_ 7-1 771 ^ I 771 771 / ' V 771 771/ 

v WW) (p'\v\p) 

2-^1 TP, _ TP, F, - F„ 



(3.24) 



E — Ei E — Ei 1 E-Ei-E(E,pY 
p'¥=p 
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W=Wp) + E Mi + .., ,,25) 

The above formula is equivalent to write the full Green's function in interacting systems 
with a self-energy part and a free part. The self-energy part is defined from the one-particle 
irreducible part in standard manner. In the present system, the one-particle irreducible part 
is defined from the momentum, and the self-energy part T,(E,p) is defined based on it. It is 
easy to see it in diagrams. They are given in Fig. 3. 

In the above perturbative expansion, if the energy denominator E — E\ does not vanish, 
there is no infra-red divergence. Moreover the series converge fast if each term in Ti(E,p) is 
small. In fact, the boundary potential V has a non-zero value only in narrow regions near 
the boundaries. Hence the matrix element from the boundary potential, (p\V\p') is of order 
Voa/L. There are the following Ai and A2 : 

For such E that satisfies 

\E-EA > A 1; (3.26) 

H(E,p) satisfies, 

Z(E,p)<A 2 , (3.27) 
A 2 < Ai, 

where Ai is of order Voa/L. The convergence of Eq.(3.25) is good then, and S(p) has no 
singularity. Hence the energy region of Eq.(3.26) is regarded as QHR. They have the width 
of order a/L and becomes very small if L is large or the magnetic field is strong. It should be 
noticed that Ai is finite and QHR exists if Ai is smaller than Landau level energy spacing 
£7+1 — Ei. This is possible even in a finite L and a finite magnetic field case. 

Sq{p) in free systems are replaced with S(p) in systems with disorders. In QHR, S(p) has 
no singularities and the perturbative expansions due to disorders which modifies So(p) into 
S(p), converge well. Their effects in the a xy are also studied by perturbative expansions and 
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they do not give any corrections to the quantized Hall conductance from the Coleman-Hill 
argument. Consequently the quantization of the a xy is exact in QHR that is defined by 
S(p) of Eq.(3.23). The fact that the QHR is defined by the momentum representation, S(p), 
becomes important when we discuss the edge states generated by the boundary potential. 

(iv) Systems with interactions in quantum Hall regime, \Ep — E[\ > 5. 

Since all the one-particle states are localized in the energy region, \E — Ei\ > 5, current 
correlation function has no singularity if the Fermi energy is in this energy region. Since 5 
and Ai should be the same order, we identify them. We investigate this energy region here 
and show that the interactions do not give any corrections to a xy in QHR. 

As was shown in the previous part, the momentum non- conserving part of the current 
correlation function has no first derivative at the origin and does not contribute to a xy . A 
momentum conserving part in a system with interactions also has no first derivative at the 
origin if that agrees to a (regular) limit of a momentum non-conserving part. 

In QHR, TT^ u (q) has no singularity in q^. Perturbative higher order diagrams due to 
interactions and disorders also have no singularity. They are computed by integrating a 
product of propagators and vertices. By cutting one of bosonic propagators and regarding 
them as two external lines with different momenta, following Coleman and Hill 12 , we are 
able to regard the original amplitude as a limit of momentum non-conserving amplitudes, 
n^v '(qi, 52), which satisfy the same identity as ^1/(91,52) of Eq.(3.14). Hence, tt^u(qi,Q2) 
has no linear term in qjj) and n^qi, q2)\q 2 ^qi does not contribute to a xy in QHR, as well. 

Only the momentum conserving part that is isolated from the momentum non-conserving 
parts and is not connected with any non-conserving parts contributes to a xy . There are two 
kind of momentum conserving parts that are isolated from non-conserving parts. The first 
type is the lowest order diagram, i.e., one loop diagram. This is obviously isolated from any 
other diagrams. The second type is similar to the first one, but is a loop diagram with a 
dressed momentum dependent propagator. The momentum dependent propagator expresses 
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extended states and has a singularity within a finite band width, which is estimated later. 
If the Fermi energy is in the outside of the continuous energy band region, all the particle 
states around the Fermi energy are localized and the system is regarded as QHR. 

Thus the higher order terms do not contribute to a xy in QHR. If there are ultra-violet 
divergences, physical observables are computed with renormalized quantities. The bare 
charge is replaced with the renormalized charge. Consequently, a xy is given by exactly 
quantized value (e 2 /h)N in QHR with the renormalized charge e. 

§4. Cylider geometry 

We study electron system of a cylinder geometry with a length L. For simplicity we 
assume L is an integer multiple of lattice spacing a. Wave functions in center coordinates 
Rmn satisfy, 

V(R m , n + Le x ) = q(R m , n ), (4.1) 

and the corresponding discrete momentum p x satisfies, 

Lp x = 2nn, n = integer. (4.2) 

If the L is much larger than magnetic lattice spacing a, localized one particle states are 
essentially the same as those of the previous infinite system, because they have finite spatial 
extensions and are insensitive to the boundary conditions. Hence one-particle states in the 
energy region \E — E{\ > 5 with a slightly modified 5 are localized. This energy region 
corresponds to QHR. We study a finite size effect of a xy in QHR. 

The Hall conductance in QHR of a cylinder geometry is given by, 

e 2 f 1 dS~ l -dS -1 -dS -1 

°*y = ^22j dp°d Py e^^- 2 Tt[ S_S S], (4.3) 

Px 

where the p x integration in Eq.(3.7) was replaced with the discrete momentum summa- 
tion. By the replacement, the integral is generally changed depending upon the form of the 
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integrand. In the present case, however, the integrand becomes, 



1 



Tr[ 




dpn dp v dp p 



4tt 2 eB 



i 1 



E 



po - Ei ± ie 



1 



(4.4) 



where — ie(+ie) is taken if the energy E\ is less than(larger than) Fermi energy, and does not 
depend on p x . Hence a xy is unchanged. a xy in the cylinder geometry is identical to that of 
the infinite system. There is no finite size correction in quantized Hall conductance in QHR. 

The effects of impurities and interactions in QHR are studied in the same way as infinite 
systems. They do not give any corrections. We will describe this in more detail in the next 
section. 



We study realistic systems which have Hall bar geometry in this section. 
(5-1) finite system 

Realistic experiments are done with electronic systems of Hall bar geometry. In one 
direction a Hall bar system has potential barriers which confine electrons in inside. There 
are gates in perpendicular direction through which electrons move in and move out. So it 
would be sufficient to study a system that is finite in one direction due to potential barrier 
and infinite in another direction. 

The potential barrier in the positive x region is expressed with the following potential 
term : 



The potential is invariant under translation in the y-direction, hence the eigenstates are the 
plane wave in this direction. The eigenvalue equations were solved before in the present 
representation. The edge states which have the p y dependent continuous energy in the 



5. Hall bar geometry 




(5.1) 



V(x) = Vo(0(-x)+6(x-L)). 
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region E\ < E < E\ + Vq and are localized in the x-direction are found. They cross the 
Fermi energy. Since edge states have the continuous energy around the Fermi energy, they 
could contribute to conductance, generally. However, they are confined along the edges 
and move in one direction without reflected. We see that they do not contribute to the 
conductance, then. 

The one-dimensional edge states near the Fermi energy are chiral and electrons move in 
one direction. Their properties are unchanged by adding disorder potentials, because in the 
wave equation for a chiral mode with a potential term, 

[i-j^ + V(x)]V = EV, (5.2) 

the potential can be removed completely by a gauge transformation 13 , 

d - ~ (5-3) 

i—m = Em. 

ox 

Hence the phase factor of chiral modes are modified by the potential but the waves are neither 
reflected nor localized. They carry electric current but do not contribute to conductance in 
the above situation. 

We use the resummation formula in the momentum representation and study the finite 
size effect. If one edge state at one side interacts with one edge state at another side, 
extended states can be formed, then the propagator in the momeutum representation has 
singularities. The effect of the edge states and the boundary potentials are found by writing 
the Hamiltonian as, 

H = ^2bi 1 (m 1 ,n 1 )ai 2 (m 2 ,n2){Ei 1 5i u i 2 5 mi ,m 2 Sn 1 ,n2 + Vi 1 j 2 (m 1 ,n 1 ;m 2 ,n 2 )}. (5.4) 

The transformed potential V/ 1; / 2 (mi, n\\ rri2, n<2) vanishes in inside region and becomes diag- 
onal form ^ ^m^m^ni,^ i n outside region. In narrow boundary region of few magnetic 
lengths, it has off diagonal term. The potential and eigenfunctions are transformed in the 
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y-direction as Eq.(2.3), and satisfy Eq.(2.4). We solved the equations numerically. Energy 
eigenvalues and corresponding eigenfunctions are given in Fig. 2. Ovbiously, edge states are 
localized in narrow region in x-direction and have energy in a range E\ < E < Ei + Vq. This 
is understandable easily from the properties of V/ l5 / 2 (mi, n\; 1712, 112)- 

Edge state at one boundary could couple with edge state at another boundary if their 
distance is small and there are impurities and interactions. They could give singularities 
in the momentum representation, S(p), then, and contribute to the conductance. To study 
their effects, it is convenient to write the previous Hamiltonian, Eq.(5.4), in the following 
manner : 

H = Hq + Hi + H2, 

X inside 

H 1= (E l + V )b l (X)a l (X), (5 ' 5) 

X outside 

H 2 = miuh(XiiX2)b h (Xi)a l2 (X2). 

boundary 

Single particle energy from Hq and Hi are E t or E\ + Vq respectively. (5V)i 1 j 2 (Xi, X2) is 
not vanishing only when Xi and X2 are in the narrow boundary regions of few magnetic 
distances, hence a perturbative treatment of H2 is possible. When the Fermi energy is 
slightly larger than E\ but is much smaller than Ei + Vq, the one particle eigenstates of Hi 
decouple and contribute to conductance only through virtual effects. Hence we take into 
account Hq and H2 meanwhile. H2 is taken perturbatively and virtual effects are studied 
later. We study particles confined in the finite inside region. 

We express the field operators of the Landau level in momentum representation as, 

fir/a 



„ J-tt/o, a 



k(X) = I a dpy-e^kif), 

J-Tv/a a 



(5.6) 

n/a 



2ir 



Px = ~j—n x + a. 

J^x 
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In Eq.(5.6), n x is an integer and a parameter a depends on boundary condition. The current 
operator, J M , which is expressed with these operators, is conserved within the Hilbert space 
of confined particles. It satisfies, 



The current operators, J^, satisfies the same commutation relation as that of infinite system, 
Eq.(3.4). Hence the identities of the infinite system, Eq.(3.5), are satisfied, but slight modifi- 
cations are necessary. Current correlation function is expanded with the momentum carried 
by current, which becomes discrete, and its slope at the origin, which is simple derivative, is 
proportional to the conductance. One loop diagram contributes in QHR, from Coleman-Hill 
theorem. The momentum has a discrete component, from Eq.(5.8), hence the Hall conduc- 
tance agrees to that of the torus geometry, Eq.(4.3), which has no finite size correction and 
is independent of L x . 

Coleman-Hill theorem in finite systems. 

In the present geometry the x-coordinate is defined in a finite region and the correspond- 
ing momentum, p x becomes discrete. Hence the derivations of Ward-Takahashi identity and 
related low energy theorems should be re-examined. 

We define the current correlation function, vertex function, and the propagator as 
Eq.(3.3). For a convenience, we use a notation j^(x) instead of j^(x) of the previous part. 
We concentrate to the cases where Green's functions in configulation space vanish in the 
outside region due to the finite Fermi energy that is smaller than Vo. Integrations of the x 
coordinate, then, are defined in a finite region, < x < L, and the corresponding momenta 
thus become discrete from completeness. 

The current expectation value is connected with external vector potential through the 



•9% = C, 




(5.7) 



X or Y in outside 



(Confined particle \C\ Confined particle) = 0. 
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current correlation function as 

J ll {x)= \ dyir IIV (x,y)A„(y), 

J , (5-8) 
7r^(pi,p2) = J dx 1 dx 2 e ip ^ +ip ^7r^(x 1 ,x 2 ). 

The current correlation function has a translational invariant term and a non-invariant term. 
They are expressed by 

^J(pi,p 2 ) = (2n) 2 S^\p 1 +p 2 )LS^ +p , n^(p 1 ), 

1 2 ' (5.9) 

4?(PbP2) = (2 7 r)5( 1 )(p?+p°)4?(p 1 ,p 2 ;p?), 

or 4?( Pl ,p 2 ) = (2tt) 2 5( 1 )( P ? +^)5W(pf +pi)4?(Pl,Pf), 
where <5(x) is Dirac delta function and 5 m> o is Kronecker delta. If the Fermi energy is in the 
energy gap region or in the localized state region where all the energy eigenstates around 
the Fermi energy are localized, (p) and nffl (pi,p 2 ) have no singularities. They satisfy, 

/45(p)=/45(p) = o, 



(5.10) 

p°A 2 J +rf4 2) = p°2^(pi,p2)+p\^ = o, 

and are expanded with momenta as 

= Ce UVD p^ + higher power, 

(5.11) 



tt$(pi) = Ce^ p p{ + higher power, 



7r$(pi,p2) = C(p 2 ), 

from the arguments of Section 3. Note that the coefficient C and higher power terms are 
defined uniquely by the slope or the higher order curvaures at the origin and are defined 
by the simple derivatives, even though the momentum is not infinitesimal quantity in finite 
systems. We see that C contributes to the Hall conductance, because the total current is 
given by, 

J dx dx' 2 £ dx[J y (x[) = I d ^^-\ X>S(°> P2) +45(0,P2)M,(p 2 ). (5.12) 

P2 

At a point of measurement there is no electric field, because that is a gate, and derivative 
of A v vanishes. Hence, ^v{p\-,p 2 ) and higher terms of nfy{p) thus do not contribute to the 
Hall conductance. 
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Higher order correction due to interactions in quantum Hall regime does not exist in C, 
as in the infinite system. Higher order corrections of 7r^(p) come from diagrams in which 
all the internal momenta satisfy boundary condition. If an internal momentum does not 
satisfy boundary condition of Eq.(5.6), the momentum in that direction is not conserved. 

(2) 

Consequently, this kind of higher order diagram contributes to only (pi,p2), and does 
not contribute to C. Examples are given in Fig. 7. 

As was done by Coleman and Hill, in order to study ^}{p) we study new diagrams in 
which internal lines of tt$(p) are cut and two different momenta, l\ and I2, are given to these 
lines. Whenever h + fo ^ 0, two current have different momenta. These diagrams then cor- 
respond to 7r^(pi,p2), and are written as ir < $(pi,p2) = 5^\pi+p2 + h + h)^^{pi,P2), and 
do not contribute to C from the above argument. Amplitude 7t^(pi,p2) have no singularity 
in QHR and linear terms in p\ or p2 are not allowed, because p\ and P2 are independent. 
Then, n^(p,p) has no linear term also. Thus, none of higher order diagrams contribute to 
C and only the one loop diagram contributes to C. 

The linear coefficient is written as a topological invariant. Since the current conservation 
and the commutation relations are the same as the infinite system, the transformed vertex 
part as Eq.(3.8) and (3.9), hence satisfies, 

q»f^p u p2) = S- 1 ^) - S-\ P1 + q). (5.13) 

Consequently, the derivative is obtained by comparing a linear coefficient of both sides 

as 

The above equation is satisfied with the simple derivative in the left hand side, despite of 
the fact that one component of the momentum is discrete. Using this formula, the linear 
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coefficient of the current correlation function is expressed as, 

.2 Q 



>xy 



3! 



dp p 



' 1 E /^CCC* 



(5.15) 



This formula gives topologically invariant expression of the Hall conductance in finite sys- 
tems. 



(5-2) Boundary effect 

We study off-diagonal term H2 of Eq.(5.5). This term does not conserve momentum in 
the x-direction. Hence perturbative expansion of this term gives higher order diagrams in 
which the total momuntum is not conserved in the x-direction in addition to the momentum 
conserving diagrams. These higher order diagrams do not contribute to the a xy as far as 
they are treated perturbatively. Perturbative treatment is good if the Fermi energy is in 
the outside of the continuous band of the momentum conserving propagator, but may not 
be good if the Fermi energy is in the inside of fictitious band of section 2. Real extended 
states could be formed by boundary potentials in a finite energy region around the center 
of Landau level. The width of the fictitious band in our momentum representation can be 
made small in certain conditions discussed before. The QHR is in the outside region of the 
energy bands of real extended states and of fictitious extended states. The width of the 
bands was estimated in Section 2 and was given by Eqs.(2.9) and (2.11). Suitably large 
L x or strong magnetic field make the widths much smaller than the Landau level energy 
spacing, E\ — E\_\ = eB/m. There are enough spacings for the localized states, then, and 
QHR do exist in this situation. 

In QHR, the momentum non- conserving term gives no contribution to the slope of cur- 
rent correlation function from the arguments of Eqs.(5.11) and (5.12). if 2 in Eq.(5.5) does 
not conserve momentum in the x-direction but conserves in the ^-direction. Consequently, 
impurities and interactions do not modify the value of the quantized Hall conductance there. 
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In the inside of extended states energy band of S(p), the expression Eq.(5.15) is valid 
but there are corrections to the quantized a xy and to a xx . Their contributions were found 
before, and are given by 



e 2 



2tt 



e ~ e, 



(5.16) 



where e and e' are small parameters that are proportional to the number of the extended 
states. 



(5-3) Virtual effects (renormalization effect) 

Higher energy states, such as higher Landau level states, and the states in the outside 
region expressed in H\, give only virtual effect in higher order corrections. In the infinite sys- 
tem, by using the Ward-Takahashi identity derived from the current conservation, the charge 
renormalization factor cancels exactly with the vector potential renormalization factor 1 and 
the final formula of the Hall conductance, 

N w , 

(5.17) 



e 

®xy = ~J~^wi 



1 ( m r as_1 (P)s, .dS-^lp)-, .dS-^p)-, „ 

was given. Now in a finite system, the current conservation and the Ward-Takahashi identity 
are satisfied. Hence we have the formula, 

°xy — ^ J^ w , 

AT' 1 f X- T r^" 1 Q 9 ^ Si (5 ' 18) 

By combining Eq.(4.4) and Eq.(5.18), we have the Hall conductance in a QHR of finite 
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system as, 



e 2 



a xy = —N, N = integer (5.19) 



The a xy has no finite size corrections in QHR. 
(5-4) Finite current effects. 

The linear response formula, Eq.(5.8), is valid actually if the current and the vector 
potential are infinitesimal, 

5J^{x) = J dyir^(x,y)5A v (y). (5.20) 

In order to know the relation for a finite current, we integrate the above relation. That 
is possible once the current correlation function is computed in a system of a finite vector 
potential. Our formula and theorem can be applied to such correlation function. Then 
the infinitesimal current, infinitesimal potential, and current correlation function depend 
inplicitly on the vector potential. We parametrize them with a parameter s, under boundary 
conditions 

Jn(x, 0) = 0, J^(x, 1) = J M (x), 

A^x, 0) = 0, Afj,(x, 1) = A^x), (5.21) 

5J^(x,s) = J dyn^(x,y; A^(s))SA v (y, s). 

A total current satisfies a similar equation, 

SI x (s) = a xy [A ll (s)\5Vy(s), 

Vy(0) = 0, Vy{\) = Vy, (5.22) 

7^(0) = 0, I x (l) = I x . 

Now, the a xy could have the same value during the change from s = to s = 1, if the 
electronic system is in the same QHR. Then the relation for a finite current is obtained by 
integrating Eq.(5.22) with s, and becomes 



e 2 



Ix = CxyVy, a xy = -j^N- (5.23) 



Note that this is possible only in the plateau region 14 . On the other hand, if the system 
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moves from one QHR to extended states region or to another QHR, then we are not sure 
if linear relation between a finite current and a finite voltage is satisfied. The relation may 
become complicated, then. We will give more arguments of finite current effects, especially 
on a breakdown of the QHE, in the next section. 

The relation between the finite total current and the finite total voltage becomes linear 
in the plateau regions, and the conductance a xy is quantized. It behaves as Fig. 8. 

§6. Current distribution, Biittiker-Landauer 
formula, and Breakdown of QHE. 

In our proof of the quantum Hall effect given in the previous parts, current distribution 
is irrelevant. The Hall conductance is the ratio between the total current in one direction 
and total voltage in another direction and so is quantized exactly in QHR. The current 
density is not uniform and varies with spatial region, generally. We study its implication in 
this section. 

There is a completely different approach of the quantum Hall effect from ours. In its 
proof, it is assumed that edge states are the only current earring states around the Fermi 
energy and is used Biittiker-Landauer formula for one-dimensional systems. The formula 
may be valid only under the assumption. The current distribution is thus important. Two 
approaches give totaly different result when the current becomes larger. We study the current 
distribution and finite current effect in this section. We will see that in general situation, the 
bulk states as well as the edge states carry the current. Hall electric field is thus generated 
at the bulk and it leads Landau level broadening. QHR becomes narrow consequently and 
vanishes eventually. Breakdown of QHE occurs. 

(6-1) Current distribution 
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In order to compute the current distribution, we compute an effective action of an 
external vector potential. The vector potential is regarded as either the external potential 
added to the system or as a Lagrange multi-plier which is expressing a condition of a finite 
external current. Variational principle to the effective action gives the current distribution. 

We study the system in which the vector potential couples with electron field in a gauge 
invariant manner as, 

J dx[tf(ih^ + eA )V - ^t (P + e ^o + e ^) g _ qr^V b (x)\, 

|V x A \ = B, (°) 

Vb : boundary potential. 

We integrate the electron fields ^f(x) and ^(x) and obtain the effective action of Aq and A. 

Long distance phenomena are represented by low dimensional terms of the effective action 15 , 

\ J dx[c(x)Fl + e(x)F$ y - d(x)F$ + °-f J dxe, vp A^ A? + • • • , 

Fq x = doA x — 8 x Aq, Fo y = doAy — d y Ao, F xy = d x A y — d y A x , 
The x-integration is defined in a finite region from to L. The coefficients c(x), e(x), and 
d(x) of the above action are the constants in the bulk and change their values near the edge. 
They are computed as, 

e 2 m 

c(x) = e(x) = — — , 

3g2 (6.3) 

in the bulk. Near the edges they change the values. Because we need only c(x) here, we give 
c(x) near the edges for a boundary potential with a constant slope, E^. 



e2 m „2 B I 2 „-2eBAx 



„2 



c(x) = — — - e z — d—e- zeti * x , (6.4) 
v ; 8vr eB AE h V eBn ' v ; 

where Ax is the distance between the boundary and the coordinate x. The first term comes 

from inter Landau levels and the second term comes from intra Landau levels. They have 

opposite sign. For high potential barrier, the first term is dominant and c(x) is positive 

definite in all the regions. We study this case here and discuss general cases in Appendix C. 
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The total current in the x-direction is given by, 



I y = f dxJ y = a xy {A°(L)-A°(0)}. (6.5) 
Jo 

We assume the translational invariance in the y-direction and obtain a stationary solution 
of the effective action under the constraint of total current I y with time independent ansatz 



d Ai = 0, d y A = 0. (6.6) 
Then we minimize the following action, 

- J dx[c(x)(d x A ) 2 ] - fia xy J dxd x A (y 2 ). (6.7) 

Euler-Lagrange equation is given by, 

d x (c(x)d x A ) = 0, (6.8) 

and is solved as 



c(x)d x Ao = const = Co, 

C (6.9) 
d x A = — -, c(x) ^ 0. 
c(x) 

The constant Co is given from the condition of total current, 

Co \ dy—— = — ^- = const. (6.10) 

Jxi c \ x ) a xy 

Obviously, if the coefficient c(x) were constant, the local electric field, as well as the local 
current density, would be uniform. On the other hand, if c(x) varies, electric field varies 
also. In fact from Eq.(6.4), c(x) decreases toward the edge and the electric field and current 
density increases toward the edge. The edge current is only a portion of the total current. 
In an exceptional case, where c(x) vanishes or becomes negative at the edge and stays at the 
constant value in the bulk, there is an energy minimizing solution which has a current only 
at the edge region. In this situation, it is obvious that Biittiker-Landauer formula could be 
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applied. Our formula is applicable to general situations. A general discussion concerning a 
connection between the sign of c(x) and the current distribution is given in Appendix C. It 
will be shown that edge current states may make a transition to bulk current states. 



(6-2) Finite current effects and breakdown of QHE 

Potential thus obtained modifies one-particle properties. Their effects become important 
if the magnitude of the total current is not infinitesimal but is finite. We study the system 
of the uniform electric field. Landau levels with the uniform electric field are not degenerate 
in energy and have a finite width, as is given in Appendix D. If the width in the momentum 
representation exceeds the Landau level energy spacing, QHR disappears. Then, breakdown 
of QHE occurs. 

In the von Neumann lattice representation it is easy to find the width in the momentum 
representation. From Eq.(D.2), the gauge invariant width is given by, 

eEa. (6.11) 

QHR disappears, when the band width exceeds the Landau level's energy spacing, Hcu c , for 
all Landau levels without spin effect or for energy splitting from Zeeman energy due to 
magnetic moment //, \iB. The critical electric field satisfies, 

eE c a = hu c , (6-12) 



or 



eE' c a = fiB, (6.13) 
= m" " = \7B = ^ (6 - 14) 



where Ib is the magnetic length used usually. Thus, the critical electric field, E c for Landau 
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level splitting and E' c for spin splitting are given by, 

E c = = N1B 3 / 2 , 
ea 

= flB = 3/2 

c ea (6.15) 
Aq = 25.4 x 10 3 [Vm _1 T" 3 / 2 ], 

N 2 = 0.8Ag x 10 3 [Vm _1 T" 3 / 2 ], 

In both cases the critical electric fields are proportional to B 3 ! 2 and the numerical constants 
N\ and N 2 are constant and independent of Landau levels, g-factor extracted from the 
experiment 19 is 7.3 and N 2 becomes comparable to N\. The critical electric fields for even 
plateaus and odd plateaus are computed from the previous values as, 

E° dd = N odd B 3j/2 , 



(6.16) 

N odd = jsj 2 = QA8 x i0 3 [Vm _1 T" 3 / 2 ], 

jyeven = ^ _ ^ = 1Q Q x 10 3 [Vm™^" 3 / 2 ] . 

These critical fields are compared with the recent experiments of Kawaji et al 8 . Both of 
them have _B 3 / 2 behavior and are independent from Landau levels. Our values of iV odd and 
7V even , however, are much larger than the experimental observations. 

The critical electric field has been estimated before from naive overlappings of wave func- 
tions by Eaves and Sheard 20 and from semi-classical method by Trugman, and Nicopoulos 
and Trugman 21 . In the former method, the critical electric field behaves as B 3 ! 2 , and its 
magnitude thus obtained was similar to the current values. In the latter method, the critical 
electric field behaves as B instead of £> 3 / 2 . We estimate the band width of the extended 
states due to Hall electric field first and compute the critical electric field from a condition 
that the bands overlapp. Our method is thus quite natural and leads to reasonable qualita- 
tive agreements, in the -B-dependence and the level dependence of the critical electric fields 
with the experiments. 
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Edge current systems make transition to bulk current systems, from Appendix C. Con- 
sequently, the systems of only the edge currents should show the breakdown of the QHE in 
two steps. In the first step, an edge current system becomes to a bulk current system, and 
in the second step, the QHE is broken. 

§7. Summary 

In the present work, we have shown that the quantum Hall regime(QHR) is realized in 
finite two-dimensional electron systems if the magnetic field is strong enough and that the 
quantized Hall conductance has no finite size correction in QHR. They are shown by the 
use of magnetic von Neumann lattice representation, which has been used by us before and 
is quite useful for studying one-particle properties and for showing the connection of Hall 
conductance with the topological invariant and the absence of corrections in quite general 
systems with disorders and interactions. The momentum representation is used, and QHR 
is defined based on the momentum representation of the propagator. 

In magnetic von Neumann lattice representation, base functions and dual base functions 
are local functions and have values around rectangular lattice coordinates. It is easy to apply 
the momentum representation in this method. Then Ward-Takahashi identity is expressed 
with simple and transparent form by use of multi-pole expansion technique. The derivation 
of the exact low energy theorem about the a xy is given based on them. From the theorem, the 
a xy at the QHR is quantized exactly as (e 2 /h)N and has no correction from finite size effect, 
disorders and interactions. The edge states are extended along the boundary and have 
continuous energies across Fermi energy. Nevertheless they have small overlapp with the 
momentum eigenstates, and the QHR is realized at the outside region of the singularities of 
the propagator in the momentum representation, S(p). Edge states carry the electromagnetic 
current together with the bulk extended states. Both states contribute to Hall conductance 
and Hall conductance is quantized in QHR. Hall conductivity, on the other hand, may depend 
on spatial region and is not quantized generally. 
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We studied the current distribution in Section 6 and found that the bulk states as well 
as the one-dimensional edge states carry the current generally. In QHR, the bulk extended 
states have energy gap but the edge states have no energy gap. Around Fermi energy, there 
are only one-dimensional edge states. Biittiker-Landauer formula may be applied to them, 
then. 

The current in QHR does not cause energy dissipation. One-dimensional chiral modes 
near Fermi energy carry a current without energy dissipation and two-dimensional extended 
states at the bulk with finite energy gap also carry a current without energy dissipation. In 
our approach, we have studied combined total current and total voltage, and we have shown 
that the a xy are quantized exactly. 

In Biittiker-Landauer approach, only the states near Fermi energy are studied and it 
was shown that their contributions to a xy gives the quantized a xy . We discussed when the 
current flows only in the edges. These systems, however, are changed to bulk current systems 
if the current exceeds the critical value. 

The current in the bulk produces the Hall electric field in the bulk. Due to the electric 
field, the bulk one particle states are modified to become extended energetically, and QHR 
becomes narrow and eventually disappears at a critical current. A theoretical estimation is 
made. The results are consistent with the recent experiment by Kawaji et al and others 8 . 
The critical electric fields are proportional to £> 3 / 2 and are independent of Landau levels 
in consistent with the experiment by Kawaji et al and others 8 , but their magnitudes are 
substantially smaller. 
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Appendix A 



The explicit form of matrices of Eqs.(3.5) and (3.6) are following: 

ad d 
dxip) = d' x (p) = -^(q^ ~ ^ga(p x ,p y ), 

o,0 

dy{p) = d' y {p) = -^dp^ + l0 g«(^'^)' 



dx{p)h,i 2 = d' x (p) hM = -jJ==(\fr 2 5 hM ^i - y/T 2 + lS bt i 2+1 ), 



1 



d y(p)h,h = d ' y (p)h,h = ~^§^ 26h ^- 1 + V 7 ^ + ib~ hh+1 ), 
di satisfy the commutation relation: 



d x (p) is given in Fig. 9. 



(A.l) 



a{p x ,p y ) =J2e iMmi - m ' l)+ipy{m - n[ \Ri\R[), (A.2) 



(A3) 



[d x (p)Jy( P )} = [d' x (pU'y(p)} = (A.4) 



Appendix B 

We study finite size correction of a topological invariant defined by a propagator, z 2 (p), 

^^/^- Tr C s C s C s| - (Bi) 

N w agrees with the integer if the momentum space is compact and the matrix space of 
Si ly i 2 (p) includes SU(2) group as a subspace. We see in the main text that the finite size 
correction does not appear in QHE. Finite temperature effect is similar to the finite size effect. 
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po integration in Eq.(B-l) is replaced with a discrete summation at finite temperature, and 
N w of QHE is given by, 

N, (finite size), 

E£i,(SKM)( cotll ^l /i _ El \ + 1 ^ 2, ( nnite temperature), ^ ^ 

(3 = — — , n : chemical potential. 
kT 

Finite size effect and finite temperature effect of N w depend on the form of the propa- 
gator. We compute these effects of ground state 16 of Dirac field, which has 

70 = m, 

(B.3) 

71 = m, 



72 = 2T 2 . 



N w is given by, 



Nu= lH J coth * (nnite size), ^ 
2 m [coth-^-, (finite temperature), 
where L is the width in x-direction. Thus the topological invariant N w of the Dirac theory 

has the both of finite size correction and finite temperature correction of exponential type 
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Appendix C 

We study a system described by one-dimensional Landau- Ginzburg type static energy 
for Ao, 

U = £ dx{^c(x)(d x A ) 2 + ^X(x)(d x A ) 4 }, (C.f) 

with a constraint, 

a xyly = / dxd x Ao = constant. (C.2) 
Jo 

In writting (C.f), We have ignored higher derivative terms such as (d 2 Ao) 2 . Hence a self- 
consistency of Ref.18 is not considered here. Owing to translational invariance in ^-direction, 

40 



we use one-dimensional form. We regard the coefficient c(x) and X(x) are known and study 
solutions of Euler-Lagrange equation under the above constraint. This gives the electric field 
and the current density. 

Depending on sign of the function c(x), solutions have totally different properties. 



1. c(x) > 0, \(x) > 

Euler-Lagrange equation from (C-f) is given by, 

d x [(d x A ){c(x) + X(x)(d x A Q ) 2 }} = 0. (C.3) 

Integrating (C-3), we have 

d x A Q {c(x) + X(x)(d x A ) 2 } = C = constant. (C.4) 

The electric field 8 x Aq is given by solving (C-4). For small 8 x Aq, the solution is approxi- 
mately given by, 

This solution has electric field in whole region and the current flows in the bulk and at the 
edges. The non-zero constant Cq is determined from the constraint (C-2). 



2. c(x) > 0, X(x) > in the bulk and c(x) < 0, X(x) > at the edges. 

Euler-Lagrange equation is the same as before and is given in (C-3). Because c(x) is 
negative near the edge regions, a new type of solution, which corresponds to Cq = 0, exists. 
If Cq = 0, we have 

c(x) + \(x)(d x A ) 2 = 0, 

(C.6) 

or d x A = 0. 

In a region where c(x) > 0, the electric field vanishes. In another region where c(x) < 0, 
either the electric field vanishes or the electric field satisfies the first equation of (C-6). They 
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are determined from the constraint of total current (C-2). This type of solutions have electric 
field only at the edge regions. The bulk has neither electric field nor electric current. Hence 
this corresponds to edge states. 

The edge current region where the first solution of (C-6) is satisfied is determined from 
the constraint (C-2), hence it depends on the total current. A relation between the width 
of the edge current region and the total current is given in Fig. 9. This type of solution 
disappears if the current exceeds a critical value which satisfies, 

[ L dxJ^4 = a xy J c . (C.7) 

Jo,c(x)<0 V 

In obtaining Fig.9, we have used, 

X(x) = \q = constant, 

c, at the bulk, (C.8) 

— c + dx 2 , near the edge, x < 



c{x) 



The solution becomes that of Co 7^ if the current exceeds the critical value, and is obtained 
by solving (C-4). It has current in the whole region. The current distributions of the solutions 
are given in Fig. 10. 

Appendix D 

In the von Neumann lattice representation, it is easy to compute energy spectrum in the 
momentum representation. Hamiltonian of a system with a uniform electric field, E, in the 
x-direction is given by, 

H — Ho + Hi, 

Ho — ^2Eibi(Ri)ai(Ri), (D 1) 

l,Ri 

Hi = eE^b h {Ri)[d xhh 5^ 2 + ^h,h am ^ r u r 2 + ^i u i 2 d x (Ri - R 2 )\a h (R 2 ), 

where the zero-momentum state is not included in the sets {bi(R)} and {ai(R)} from the 
constraint of the minimum coherent states. d x \ x \ 2 is given in Eq.(A.3) and the Fourier 
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transform of di(R) is given in Eqs.(3.6) and (A.l). The first term of Hi gives an inter 
Landau level mixings and the second term gives the static energy due to the electric field. 
These two terms do not contribute to the width of the Landau levels in the momentum 
representation. The last term, on the other hand, gives the intrinsic momentum dependent 
energy, d x {p). In the momentum representation, d x {p) and d y {p) are expressed as, 

a 2 d 
d x {p) = eE— Py + —6(p), 

Z7T ° Px (D.2) 

dy (p) = w m, 

with a suitable gauge function 9{p) in the momentum space. A gauge invariant quantity, 
§di{p)dp l , is used for defining the gauge invariant energy width. The invariant width is 
given by, eE^(^-) = eEa. Hence the Landau levels in systems of the constant electric field 
in the x-direction have the width, eEa. Hence the width of the extended Landau levels is 
inversly proportional to \[B. QHR disappears if the above width agrees to the Landau level's 
energy spacing, which is proportional to B. The critical electric field of the breakdown of 
QHE is proportional to B 3 ^ 2 and the proportional constant is independent of Landau levels, 
in agrement with the experiments. 
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FIGURE CAPTIONS 



1) Energy eigenvalues (a), eigenfunctions (c) at p y = n/2a, and their currents (b) are 
shown for a narrow potential well of a width, a. Left moving modes and right moving 
modes are not separated. In (c), potential barrier is denoted by shadow region. 

2) Energy eigenvalues (a), eigenfunctions (c) at p y = n/2a, and their currents (b) are 
shown for a wide potential well of a width, 5a. Left moving modes and right moving 
modes are separated. Edge states are confined in narrow edge regions. In (c), potential 
barrier is denoted by shadow region. 

3) Lowest order self-energy diagram which removes the degeneracy of Landau levels due 
to potentials is shown. 

4) Lowest order self-energy diagram which removes the degeneracy of Landau levels due 
to interactions is shown. 

5) One-particle spectrum in the presence of the edge region is shown. From Eqs.(2.9), 
(2.11), and (2.14), in narrow regions around the center of Landau levels, there are ex- 
tended states. At outside of this regions there are localized states and one-dimensional 
edge states. These are regarded as QHR. 

6) Ep dependence of the topological invariant, Eq.(3.7), in free system is shown. 

7) An example of higher order diagrams is shown. If the internal line does not satisfy 
the boundary condition, these diagrams do not contribute to (p) but contribute to 

45(Pl,P2)- 

8) Hall conductance a xy in realistic system is given. In quantum Hall regimes where 
there are only localized states and one-dimensionaly extended states, a xy agrees with 
(e 2 /h)N. 

9) The width of edge current solution is obtained from the constraint Eq.(C-2), for the 
second case where c(x) becomes negative toward the edges. Edge current solutions 
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disappear if the current exceeds the critical value which satisfies (C.7). 

10) A function c(x) of (C-8), and the current distribution for a small current case and for 
a large current case are shown. The current is restricted in narrow regions near the 
edges if J y < J c but the current spreads into the whole region if J y > J c . 



48 



